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Abstract 

We consider a simple model for a superlattice composed of a thin magnetic 
film placed between two bulk superconductors. The magnetic film is modelled 
by a planar but otherwise arbitrary distribution of magnetic dipoles and the 
superconductors are treated in the London approximation. Due to the linear- 
ity of the problem, we are able to compute the magnetic energy of the film in 
the presence of the superconductors. We show that for wavevectors which are 
much larger than the inverse London penetration depth, the magnetic energy 
is unchanged with respect to the film in free space, whereas in the case of 
small wavenumbers compared to the inverse London penetration depth, the 
magnetic energy resembles the energy of a distribution of magnetisation in 
a two dimensional space. Possible experimental applications of these results 

are discussed. 
75.70.-i, 74.80.-g 
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I. INTRODUCTION 



The interplay between superconductivity and ferromagnetism in bulk materials has been 
the subject of active research since 1957 when GinzburgS published a paper in which he 
considered the effect of the field created by a bulk distribution of magnetisation on a su- 
perconductor, which was described by the London equations!. He has concluded that for 
a ferromagnetic induction field of the sample larger than its superconducting critical field, 
this field would destroy superconductivity, but he also pointed out that in thin films or wires 
where the induction field is much smaller (due to demagnetisation effects) and the critical 
field higher (due to the small diamagnetic energy) than in bulk superconductors, it should 
be possible to observe the coexistence of the two phenomena. Experiments carried out by 
Mathias et al.i on Lanthanum with several rare-earth paramagnetic impurities dissolved at 
low concentrations suggested that the interaction responsible for the depletion of the su- 
perconducting critical temperature of Lanthanum is the exchange interaction between the 
paramagnetic impurity spins and the superconducting electrons. This interaction induces 
an effective ferromagnetic interaction between the (antiparallel) spins in the Cooper pair, 
which tends to destroy it and hence destroy superconductivity. Anderson and Suhli have 
shown that the RKKY interaction between the ferromagnetic spins due to the conduction 
electrons is significantly reduced in the superconducting state, but pointed out that ferro- 
magnetism could coexist with superconductivity if the ferromagnetic atoms formed small 
domains. The dependence of the superconducting critical temperature on the concentra- 
tion of magnetic impurities due to exchange scattering of electrons from these impurities 
was addressed with the microscopic theory of superconductivity by Abrikosov and Gor'kovi. 
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de Gennes and Sarmau have estimated that tipically, the exchange interaction between lo- 
calised moments and superconducting electrons would be 10 2 — 10 3 larger than the dipolar 
interaction considered by Ginzburg. The detailed form of the Landau-Ginzburg theory of 
ferromagnetic superconductors was worked out by Suhfl. Despite their conflicting character, 
superconductivity and ferromagnetism are seen to coexist in bulk systems, e.g. in HoMo 6 S 8 
and ErRh 4 B 4 i. 

Another system in which the coexistence of superconductivity and ferromagnetism has 
been observed in the bulk is the nuclear magnet Auln 2 This compound shows a super- 
conducting phase transition at T c = 207mK and an ordering transition to a ferromagnetic 
state at an even lower temperature T M = 35/i-fT. This particularly low temperature can be 
explained by the weakness of the interaction between the nuclear spins (which is primarily 
due to an indirect exchange via the conduction electrons). 

On a different perspective, the development of epitaxial growth of crystals has permitted 
the creation of artificial superlattices composed of superconducting and ferromagnetic ma- 
terials, e.g. Fe/V, Ni/V, Ni/Mo, EuS/Pb, EuO/Al and Nb/Gd&H. In these superlattices, 
one can experimentally study the interaction between superconductivity and ferromagnetism 
when these two effects occur in neighbouring spatial regions and also study the supression of 
superconductivity as a function of the relative proportion (i.e. layer tickness) of the two ma- 
terials. More recently, the cuprates RuS^GdC^Os-J^Kl and RuS^Gdi+^Cei^C^OioEil 
have been found to show superconductivity and a ferromagnetism below their critical tem- 
peratures, T c = 15-40-K" for RuSr 2 GdCu 2 8 _ ( 5 and T c = 37K (for an optimal x = 0.2) for 
RuSr 2 Gdi +x Cei_ 1 .Cu20io, the Curie temperatures for magnetic ordering being Tm ~ 133K 
for the first compound and Tm ~ 100K for the second. The experimental analysis shows 
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that these materials, like all cuprates, have a layered structure and that superconductivity 
and ferromagnetism seem to occur in different layers. However, a detailed analysis has been 
hindered by difficulties with the growth of single crystals^. 

Motivated by such experiments, in which magnetism and superconductivity are seen to 
occur in different spatial regions of the studied materials, we wish to address the problem 
of a thin ferromagnetic layer, placed between two bulk superconducting layers (see fig. 
[I]), in which the thickness of the superconducting layers is much larger than the London 
penetration depth of the superconducting material and the thickness of the ferromagnetic 
layer is very small compared to this quantity, which is a condition that can be easily obtained 
with the modern techniques of epitaxial growthEl. In this limit, the results obtained can 
also be applied to superlattices of the two materials, given that the ferromagnetic layers are 
decoupled from one another. 

We consider in this paper a simple model system composed of a very thin ferromagnetic 
film, with an arbitrary distribution of magnetisation in the plane of the film, which is 
placed in a spatial gap of size 2d between two semi-infinite superconductors described by 
the London equations. The film is coupled to the superconductors by the electromagnetic 
interaction, i.e. we neglect the proximity effect!! and we consider the Josephson currentEZI 
flowing between the two superconductors to be zero (the limitations of these approximations 
will be discussed in section |TV|) . Having made the approximations indicated above, we are 
able to solve the problem exactly, by firstly considering the simpler problem of a single 
dipole in the spatial gap and then superimposing the different solutions, due to the linearity 
of the London equations. One can then compute the dipolar energy of the distribution 
of magnetisation. It turns out that for wave-vectors much larger than the inverse London 



penetration depth the form of the dipolar energy in momentum space is unchanged by the 
presence of the superconductors. On the other hand, for wave vectors much smaller than 
the inverse London penetration depth, the dipolar energy in momentum space resembles the 
energy of a distribution of dipoles in a two dimensional space. This behaviour can be traced 
to the Meissner effect which confines the magnetic flux lines within the spatial gap. 

One can think of several possible ways to detect the effects of this change of behaviour 
of the dipolar interaction at low wave vectors. If one were able to choose the materials 
composing a layered geometry of superconductor /ferromagnetic film/superconductor in such 
a way that the Curie transition temperature of the magnetic film to the ferromagnetic state 
is lower than the critical temperature T c of the superconductor, one should be able to 
measure the critical properties of the system at the paramagnetic-ferromagnetic transition, 
in particular such quantities as the specific heat and the magnetic susceptibility, with the 
superconductors already displaying the Meissner effect and therefore with a modified form 
of the dipolar energy. 

It was shown by Aharony and FisheiH that in a d- dimensional system, a ^-dimensional 
dipolar interaction (such as the one occuring in a bulk tridimensional ferromagnet or in the 
layered geometry superconductor / ferromagnet / superconductor) is a relevant interaction (in 
the sense of RG) near a paramagnetic-ferromagnetic transition, leading to a crossover be- 
tween the critical exponents of the short-range ferromagnet and the critical exponents of the 
dipolar system when one approaches the critical temperaturelH^l. Pelcovits and HalperinU 
have also shown that in the case of a <i-dimensional systems with a d+ 1- dimensional dipolar 
interaction (such as the one ocurring in a thin magnetic film film in free space) the univer- 
sality class of the dipolar system is the same as above. This is due to the fact that, at the 
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fixed point, the 'effective' (renormalised) dipolar coupling constant is infinite, making the 
susceptibility of the system independent of the longitudinal degree of freedom of the mag- 
netisation, which is the one sensitive to the nature of the dipolar interaction. However, in 
real systems, measurements are not taken exactly at the critical point and one always probes 
the crossover region. In this region, the dipolar coupling constant is finite and one should 
be able to detect the distinct character of the transition if the ferromagnetic film is included 
in a layered geometry with superconductors or if the film is grown in a non-superconducting 
substrate, due to the different character of the dipolar interaction at small wavevectors. 
The ideal experiment to detect such a distinction would presumably be a measurement of 
the longitudinal susceptibility using polarised neutrons&"@. Experiments done with films 
of EuS/SrS grown on a Si substrate have shownB that the low Curie temperature of EuS 
(16.5 K) is further reduced in these geometries. The authors of Ref.0 have also performed 
experiments with films of EuS/Pb, probing the transition between the superconducting state 
and the normal state in the Pb layer as a function of the applied magnetic field. Therefore, 
EuS stands as a good candidate for a material to be used in the ferromagnetic layer. It has 
the further advantage of being an insulator (see below). 

Another possibility would be the study of spin-spin correlation functions in a magnetic 
film in the ordered phase and outside the critical region. Indeed, Kashuba§ has shown that 
the static correlation functions of an XY model with 2d dipolar interactions would display a 
behaviour analogous to that of the dynamic correlations functions of the stochastic process 
described the KPZ equation in 1+1 dimensions@, for which the form of these correlation 
functions is known. An adequate experiment to probe these correlation functions at low 
momentum compared with the inverse London penetration depth (where such length is 
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typically of the order of a thousand Angstroms) would presumably be low-angle neutron 
scattering from the magnetic fluctuations in a layered geometry. Other possible experiments 
which could probe the magnetic properties of the system in the ordered phase would be the 
use of the magnetooptical Kerr effect or of the Faraday effect on samples with a single 
magnetic layer to image such a layer. 

The structure of this paper is as follows: in section U we define our model in terms of 
the geometry of the system and the equations which describe it. We also describe the type 



of boundary conditions we have to consider. In section |HI| , we present the solution of the 
equations for a single dipole and construct the solution for a general in-plane distribution 
of magnetisation by linear superposition. In section [IV], we compute the dipolar energy of 
the distribution of magnetisation and discuss the physical limitations of the model we have 
considered. Finally, in section M, we present our conclusions. 



II. GEOMETRY OF THE MODEL AND RELEVANT EQUATIONS 

The geometry of the model is as follows: an infinite distribution of in-plane magnetisa- 
tion is placed in the plane z — 0. This distribution is constituted by single magnetic dipoles, 
placed in an arbitrary fashion with respect to one another (see fig. |2]). The in-plane con- 
straint implies that all the dipoles point in a direction within the plane. Above and below the 
distribution are two bulk superconductors which extend from z = d (respectively z = —d) 
to z = oo (respectively z = — oo). The spatial gap with size 2c? is supposed to be filled with 
an insulator with magnetic permitivity fi^. The two superconductors are identical and have 
a magnetic permitivity fi (i.e. they are paramagnetic, with relative permitivity \i T = fi/fio)- 
These supercondutors are described by the London equations (see below), which imply a 
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linear relation between the current and the magnetic field. 

This linear relation allows us to consider instead a simpler problem, the one of a single 
magnetic dipole, placed at the origin of the coordinate system and oriented along the x 
axis. Once this problem has been solved, one can construct the solution for the general 
case simply by using translational and rotational invariance in the plane and by adding the 
different solutions. The linearity of the equations will guaranty that the linear combination 
is also a solution. Furthermore, a uniqueness theorem proved by London! guaranties that 
this solution is unique. 

In the spatial gap, the system is described by the following equations^!, 

Vxh = (1) 

V ■ b = (2) 

b = /i (h + m), (3) 

where the first two equations are the Maxwell equations for the magnetic field h and the 
magnetic induction b and the third equation is the constitutive relation between the two. 
For the case of a single dipole oriented along the x axis, the magnetisation m(r) = m x <5 3 (r), 
where m is the magnitude of the magnetic dipole. 

The superconductors are described by the equations 

V x h = j (4) 

b = -AV x f (5) 
b = (6) 

where the first equation is the Maxwell equation which relates the magnetic field with the 
'free' current, the second relation is the second London equation! which relates the super- 



current j s with the magnetic induction and the third equation is the constitutive relation 
between the magnetic induction and the magnetic field. The constant A is dependent on 
the type of the superconductor. In a static situation such as the one we are considering, 
the electric field e = in the superconductor and the total current j = j s , i.e. there is no 
normal component of the current. 

Substituting equations (§), @ in equation (§J), one obtains 

V 2 J - Xfj = (7) 
V 2 b - Xfb = (8) 

where = (A/^) 1 / 2 is the London penetration depth and where equation (|j) follows from 
taking the curl of (0) and using ([5]) and where we have used the fact that V • j = (equation 
of continuity) and V ■ b = 0. These two equations show that the magnetic flux density 
and the supercurrent penetrate a layer of thickness Xl at the surface of the superconductor 
(Meissner effect). 

These equations have to be suplemented by boundary conditions at the surface of the 
superconductors. These conditions are the continuity of the normal component of b, of the 
tangential components of h and of the normal component of the current j at the boundary 
surfaces of the two superconductors!. If one chooses j = V x (gz), where g(r) is a solution 
of 

V 2 £ - X L 2 g = (9) 

then one can satisfy the equations (0) and ([8]) and the boundary condition j z = at z = ±d. 
Notice that this choice implies that j z = throughout the material which is physically 
reasonable, since j z = at the surfaces z = ±d and also for z = ±oo. 
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On the other hand, in the spatial gap, we obtain from equation (|J) h = — V$m- Substi- 
tuting this result in equation (|) and using equation (0), we obtain 

V 2 $ M = V- (mxd 3 (r)) (10) 

which is Poisson's equation. Since we know the solution of this equation in free space (i.e. 
in the absence of the superconductors) to obtain the solution in this case, we can write 

, , , mp cos (j) / , x /-, -, s 

M(r) = 4tt( P 2 + + X(P ' ^ Z) ( ' 

where the first term on the rhs is the solution in free space and the function \{pi 0> z ) is a 
solution of the Laplace equation, V 2 x(p, <f>,z) =0 and where we have used cylindrical polar 
coordinates for later convenience. 

Therefore we need to solve the modified Helmholtz equation (Q) in the superconductors 
and the Laplace equation for \ in the gap and then fit the two solutions using the continuity 
conditions for b and h at the boundary. 

We can further simplify the problem if we notice that the system is invariant under a ir 
rotation around the x axis. This invariance imposes the conditions 

g-{pA,z) = -g+{p,-(f>,-z) (12) 

X(p,<f>,z) = x(p,~<f>,~z) (13) 

where g + (respectively gJ) is the solution of the Helmholtz equation in the upper (respec- 
tively lower) superconductor. Since the magnetic flux in the superconductor is given in 
terms of g by 
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and we have, in the spatial gap b = p (— V$j\/(r)+m x <5 3 (r)), then the continuity conditions 
for b and h imply that at z = d, 
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(17) 
(18) 
(19) 



where we have used the fact that g + is a solution of the Helmholtz equation and that, in 
the superconductor, b = ph. A similar set of conditions is valid at z = —d but they are 



trivially related to these conditions by equations Q12D and ([13]). 



The above equations and boundary conditions are sufficient to determine the solution of 
the problem within the London approximation. 

III. THE SINGLE DIPOLE SOLUTION AND THE GENERAL SOLUTION FOR 
AN ARBITRARY DISTRIBUTION OF MAGNETISATION 

We concluded in the previous section that in order to find the field and current distri- 
butions for the case of a single dipole, one needs to find a joint solution of the Laplace and 
Helmholtz equations, which satisfies the appropriate boundary conditions 0l7|) -(|l9|). Such a 
solution can be most easily found using cylindrical polar coordinates and is given in terms 
of Fourier-Bessel integral transforms by 
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771 i p COS (j) f°° 

$m(p, (f>,z) = — ^ + /2 + y dkk J^k p)cosh(kz) cos i 



e -kd( ,,-1 



l (^Jk 2 + Xf-k) 

x V / | (20) 

cosh (fed) + p r l J k 2 + X L 2 smh(kd) 



for $m(i*) and 



At: fi r XI J 

x 1 / (21) 

/c cosh(A;d) + /i r 1 J k 2 + A L 2 sinh(/cd) 

for g+(r) with <?-(p, 0, 2) = — g+(p, —(ft, —z) and where Ji(x) is the Bessel function of order 
one. These integrals can only be calculated explicitly in the particular case d — > 0, /i r = 1. 
We obtain 



m f 1 e~ p/x 
4tt \pA L p : 
m cos 
4^ 



$ M (p, 0, 0) = — ( — + — ^- ) cos (22) 
3+^,0,0) = ——^. (23) 



The magnetic potential $m(p, 0, 0) corresponds to the magnetic potential of a dipole which 
produces a 3d field at short distances and that at distances ^> A^ produces a 2d field, i.e. 
the field produced by a dipole in a two dimensional space. This behaviour can be traced 
to the Meissner effect, which confines the flux lines to the spatial gap and to a region of 
size Xl in each of the superconductors. Although the d — case is somewhat unphysical 
(the superconductor would have to withstand an infinite field), we shall see that as long as 
d <C Xl, this type of behaviour is essentially unchanged. 

Now, in order to generalise this solution to the case of an arbitrary distribution of mag- 
netisation, we represent the magnetisation per unit of area in the form 

m(p) = ]Tm^ 2 (p-p,) (24) 

i 
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where is a dipole situated at p { = (xi,yi) and p = (x,y). In real systems, p i will corre- 
spond to the sites of a two dimensional lattice where the dipoles are situated. The solution 
corresponding to this distribution of magnetisation is given by the linear superposition of 
the solutions corresponding to each rrij, i.e. 



-47^-^1^2)3/2 

+ y m * ; (p - Pi ] r dkkjAk \p-Pi i) cosh(^) 

*f 4tt I p - Pi | J 

e- fcd (AV 1 \A 2 + A £ 2 - /c ) 

x 



k cosh(kd) + fi r 1 \ J k 2 + X L smh(kd) 



for $ w fr) and 



g + ( P , M = ~£ 7 "Vj? A) | r dfc * I P - ft I) e-v^W(-^) 
^ 4vr/i r A L | p - pj | Jo 

x 1 (26) 

cosh(£;<i) + /i" 1 yk 2 + X L 2 smh(kd) 

for (?+(r). The function g_(r) is constructed from the single dipole solution in an analogous 



manner. Note that one cannot use the equations (|T2"D and (|T3|) because we no longer have 
the 7r rotation symmetry around the x axisB Using the magnetic potential we can 

now compute the dipolar energy of the system. 



IV. THE MAGNETIC ENERGY OF THE SYSTEM 

The dipolar energy of the system can be obtained by substituting the dipole distribution 
by an equivalent distribution of loops of current, i.e. one which will produce the same field 
distribution. The easiest way to compute the energy necessary for the formation of such a 
current distribution is to compute it with the currents of the individual loops which form the 
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distribution kept constant. The dipolar energy of the magnetisation distribution is equal to 
the energy necessary to form the current distribution when the fluxes in each loop are kept 
constant, which is minus the energy computed with constant currents and is given by@S 

Em = —Ho J d 2 p m(p) ■ h(p, 0, z = 0) . (27) 

Substituting h = — V$m> with $M( r ) given by fl2"5[) and approximating the discrete sums 
on the lattice by integrals, we obtain 



/'o 



d 2 k 2tt 



£ =sy BZ (2^T< k ' m < k >>< k ' m <- k >> x < 28 > 

(H~ 1 Jk 2 + \~ L 2 + k)e kd + (n~ 1 Jk 2 + \~ L 2 - k)e~ kd 



(/i-yfc 2 + \~ L 2 + k)e kd - (/i-yp + \l 2 - k)e~ kd 
where m(k) = £\ nij e~ tk ' Pi and where the integrals over k are over the first Brillouin zone 
of the reciprocal lattice. The expression ( p8|) was obtained in the limit in which one can 
disregard the lattice structure of the dipole distribution. If such a structure has to be taken 
into account, then one has to use Ewald summation methodsll to handle the discrete sums 
in equation (|25|). 

The expression ( p8|) has two important limits. The first is when [i r = 1 and — > oo, or 
when d — ► oo. We obtain 

8tt J bz (2tt) 2 k 

which is the familiar result for the dipolar energy of a thin film. The second limit is when 
the largest contribution to the energy comes from modes m(k) for which k < L^ 1 where L 
is a length such that L ^> ^> d (we take Al d). In this case, we obtain 
Ho f d 2 k (k • m(k))(k • m(— k)) / 1. 2 , 



E ~T^lL L ^ * " [i + - 2 W + 0M ] (30) 
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which shows that for k <C A^ 1 the energy of the system has the same form as the energy of a 
system of dipoles in a two dimensional space. Also, comparing equation (|28|) (with fi r = 1) 
with (|29|), one can conclude that the dipolar interaction has effectively been enhanced with 
respect to the simple film situation, since the fraction in (^Sj) is always larger than 1 when 
fi r — 1. This result can be easily understood from the fact that the energy given in fl2B|) also 
includes the kinetic energy of the supercurrentB 



The modified dipolar kernel F(k) = j , ■■ , = — given in fl28|) with 



k (^k 2 +\2 2 +k)e kd -(y/k 2 +\l 2 -k)e- kd 

fi r = 1, is plotted against the kernels 1/k and l/k 2 , which appear respectively, in (|29D and 



@,in figure |. 

The question now arises if one can indeed detect such change of behaviour in the dipo- 
lar interaction in artificial superlattices of superconductors and ferromagnetic materials or 
even in naturally layered systems like RuSr 2 GdCu 2 8 -«5 or RuSr 2 Gdi +:r Ce 1 _ : j;Cu 2 Oio. The 
detection of such a change of behaviour relies on the possibility of finding an experimental 
system which obeys a series of constraints. 

Firstly, the London approximation, which was used and which postulates a local relation 
between the current and the magnetic induction, is only valid sufficiently close to the critical 
temperature at which the superconducting phase transition occurs. But we believe that the 
qualitative features of this solution should be valid even when the London approximation is 
not (i.e. at low temperatures compared to the critical temperature). Physically, one should 
expect this type of behaviour as long as the superconductor displays the Meissner effect. 

Secondly, we have considered a ferromagnetic layer with zero thickness. The solution ob- 
tained can only be valid if one can neglect the proximity effect in the ferromagnetic layer, i.e. 
the induction of superconductivity in the ferromagnetic material by the superconductor^. 
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Otherwise, in a layer with finite thickness, an effective exchange interaction will be induced 
between the magnetic moments and this interaction also has to be taken into account. How- 
ever, this effect is negligible in insulators. Therefore, an appropriate material to choose for 
the ferromagnetic layer would be a ferromagnetic insulator^. We have also neglected the 
supression of the superconducting order parameter at the boundary between the two media, 
which always occurs in the presence of a film of finite thickness. In the case of a film made of 
a ferromagnetic insulator, the de Gennes boundary conditions tell us that the order param- 
eter is effectively quenched to zero at the boundary^. Since the London penetration depth 
depends on this parameter, this will mean that the effective London penetration depth will 
be greater than the penetration depth measured in the absence of the ferromagnetic film. 
Nevertheless, the London approximation remains valid (provided that we are working in a 
weak field situation, see below). The opposite limit, in which the coupling between the su- 
perconducting and the ferromagnetic layers is primarily due to the proximity effect has been 
considered, using a microscopic approach by Radovic et al@ and by Schinz and SchwablS 
using a phenomenological description. In both cases, the problem treated reduces to the one 
of decoupled thin superconducting layers embedded in a ferromagnet. Experiments done in 
superlattices of Fe/vS and Nb/Gdlll have confirmed these results. 

Thirdly, one cannot have any Josephson currents flowing through the insulating junction, 
otherwise the boundary condition j z = at z = ±d is not valid. This implies that the 
phase difference between the superconducting order parameters of the two superconductors 
is zero. This is the case for most superconducting systems in equilibrium but in certain 
junctions containing ferromagnetic materials a non-zero difference between the phases of 
the superconducting order parameter in each side of the junction, leading to a current in 
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equilibrium, has been predicted^. 

Finally, we have not considered the thermodynamical stability of the system we are 
working with. This implies that we are working in a weak field situation, i.e. the field 
produced by the distribution of magnetisation is much lower than the critical field at which 
the superconducting system undergoes a transition to the intermediate state (in the case 
of type I superconductors) or to the mixed state (in the case of type II superconductors). 
One way to achieve this condition is to choose either a magnetic material with a low sat- 
uration magnetisation or a superconductor with a high critical field (in the case of type I 
superconductor), e.g. V, or a high field H c % (in the case of a type II superconductor), e.g. 
Nb. 

In summary, despite the restrictions pointed above, we think that there is room for 
believing that one could manufacture systems composed of alternating magnetic and super- 
conducting layers where the above effect could be detected using the methods discussed in 
the introduction or others. 

V. CONCLUSIONS 

We have computed the magnetic field distribution of an arbitrary distribution of planar 
magnetisation in a spatial gap between two superconductors. The purpose of this calcula- 
tion is to provide a simple model for a superlattice of ferromagnetic materials and supercon- 
ductors where the dipolar interaction between the magnetic moments in the ferromagnetic 
system is taken into account. 

We have also computed the dipolar energy of such a system and we have shown that for 
low momenta compared to the London inverse penetration depth, the system has a dipolar 
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energy which resembles the energy of a magnetisation distribution in 2 dimensions and hence 
it leads to an enhancement of the dipolar energy compared to that of a simple ferromagnetic 
film. 

As to possible future directions of research one can point out the possibility to consider a 
treatment of the same problem in the line of what was done in referenceil, i.e. consider the 
Landau- Ginzburg theory of the layered system including the free energy of the supercon- 
ducting system and of the ferromagnetic layer and the interaction between the two via the 
electromagnetic interaction and the proximity effect, with the appropriate boundary condi- 
tions. However, in this case one has to take into account the non-linearity of the problem 
and the three dimensional character of a solution involving a non-uniform distribution of 
magnetisation in the plane. 

A slightly different model which can also be treated exactly with the methods developed 
in this paper is the one of a thin ferromagnetic layer between two bulk type II superconduc- 
tors above H&, i.e. with flux penetration in the form of an Abrikosov vortex lattice oriented 
along the z axis. In this case, the boundary conditions (0"H1) would be unchanged, but the 
equation determining g(r), equation would be replaced by an inhomogeneous equation 
where the role of sources is played by the vortices. Also, the ir rotational invariance of the 
system along the x axis is no longer present, but one can still draw useful conclusions from 
the reflexion properties of the system in the xy plane (se What makes the model much 
more difficult to solve is the fact that one cannot exclude the presence of Neumann functions 
in the Fourier- Bessel integrals which determine g(r) and $^(r), and also that in this case 
one cannot use the equation fl2"T|) as the expression for the total magnetic energy of the 
system (compare with@). One should nevertheless mention that a method for computing 
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the magnetic energy in the presence of vortices has been considered by Erdin et al.E-l for the 
case of a thin magnetic film interacting with a thin superconducting film. 
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Figure Captions 

Figure 1. Thin ferromagnetic film between bulk superconductors (schematic). 

Figure 2. In-plane distribution of magnetisation between bulk superconductors 

(schematic). 

Figure 3. The modified dipolar kernel F(k) (the continous grey plot) is plotted against 
the kernels 1/k (the long-dashed plot) and 1/k 2 (the short-dashed plot). We have taken 
H r = 1 and d — 0.1, = 1 in F(k) (in arbitrary units). It is seen that the function F(k) 
interpolates between 1/k 2 and 1/k, the crossover length being the London penetration depth 
Xl = I- 
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